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The string-net approach by Levin and Wen and the local unitary transformation approach by 
Chen, Gu and Wen provided ways to systematically label non-chiral topological orders in 2D. In 
those approaches, different topologically ordered many-body wave functions were characterized by 
different fixed-point tensors. Though extremely powerful, the resulting fixed-point tensors were 
mathematical abstractions and thus lacked a physical interpretation. As a result it was hard to 
judge if two different fixed-point tensors described the same quantum phase or not. We want to 
improve that approach by giving a more physical description of the topological orders. We find 
that the non-Abelian Berry's phases, T- and 5'-matrices, of the topological protected degenerate 
ground states on a torus give rise to a more physical description of topological orders. Using 
the Verlinde conjecture, we can even choose the canonical basis for the T- and 5'-matrices. It is 
conjectured that the T and 5'-matrices form a complete and one-to-one characterization of non-chiral 
topological orders and can replace the fixed-point tensor description to give us a more physical label 
for topological orders. As a result, all the topological properties can be obtained from the T- and 
5'-matrices, such as number of quasiparticle types (from the dimension of T or iS), the quasiparticle 
statistics (from the diagonal elements of T), the quantum dimensions of quasiparticles (from the 
first row of 5*), etc . 
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I. INTRODUCTION 

One of the most important questions in condensed 
matter physics is the description and the classification of 
different phases of matter. Landau's symmetry-breaking 
theory^ provides a very elegant answer to this ques- 
tion: different phases are characterized by different bro- 
ken symmetries. Thus by classifying all different bro- 
ken symmetries, we can get a classification of phases. 



Though extremely powerful, Landau's theory fails to ex- 
plain many new states found in experiments including 
the Fractional Quantum Hall(FQH) states. ^'"^ These new 
states found in experiments thus bring up once again the 
old question of how to classify different phases of matter. 

A closer look at Landau's symmetry-breaking the- 
ory reveals that the theory only describes direct-product 
states up to some small local perturbations (see Ref. 5). 
Since these small perturbations can only modify the 
direct-product states locally, Landau's theory can only 
describe states with "short-range entanglements". Intu- 
itively, these states can only account for a small fraction 
of all possible quantum many-body states. Thus what 
seems to be missing in Landau's theory is a large class 
of states with "long-range entanglement". Those states 
were named topologically ordered states in 1989 before 
the concept of entanglements becomes popular. ^'^ Dif- 
ferent patterns of long-range-entanglements/topological- 
orders correspond to different quantum phases. 

But what are these patterns of long-range entangle- 
ment? What are topological orders? It is hard to an- 
swer these questions since those patterns/orders are new 
concepts that do not even have a label. So to study 
topological-order /long-range-entanglement, we first need 
to invent labels or mathematical symbols for them. The 
first label/symbol invented was the ground state degen- 
eracy D(iQg{g) on torus (with genus g = 1)-^^ and other 
Riemann surfaces. ^'^ But it was clear from the begin- 
ning that this was not a very good label, since the same 
Ddegig) can correspond to many different topological or- 
ders. To obtain a better label, in Ref. 8,9, it was proposed 
to use the non-Abelian geometric phases^ ^ of degenerate 
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ground states on Riemann surfaces to characterize differ- 
ent topological orders. 

It was conjectured^ that the non-Abelian geometric 
phases (the non-Abelian part) form a complete and one- 
to-one description of topological orders. On torus, they 
are described by T- and iS-matrices, which generate a 
projective representation of the modular transformation 
on the torus. We can view T- and ^'-matrices as some 
kind of "non-local order parameter". Therefore, the con- 
cept of topological order is defined through the physi- 
cal properties of robust ground state degeneracy (called 
topological degeneracy)^'^ and robust geometric phases 
induced by the modular transformation of the degener- 
ate ground states, ^'^ just like the concept superconduct- 
ing order was introduced through physical properties of 
zero resistance and Meissner effect. This led to the es- 
tablishment of the concept of topological order (and long- 
range entanglements). Also mathematically, the T- and 
S'-matrices determine the whole theory for abelian mod- 
ular tensor categories. It is thus very desirable to develop 
a comprehensive theory of topological orders based on T- 
and ^'-matrices. 

However, the current comprehensive theory of topolog- 
ical orders is not based on T- and ^'-matrices, but rather, 
based on the string- net approach of Ref. 12 and local 
unitary transformation approach of Ref. 5. They provide 
a systematic description of topological orders (or more 
precisely, non-chiral topological orders^^) by systemati- 
cally labeling the corresponding "long-range entangled" 
many-body wave functions through a set of fixed-point 
tensors. These fixed-point tensors form a mathematical 
structure called tensor category theory. Since the many- 
body wave functions described by the fixed-point tensors 
are explicitly known, we can construct exactly soluble 
Hamiltonians to realize each topological order described 
by fixed-point tensor. 

Although the fixed-point tensors approach {ie the 
string-net or the local-unitary-transformation approach) 
is closely related to the many-body wave function, it 
is not clear if different fixed-point tensors always corre- 
spond to different quantum phases. So to better identify 
the quantum phases described by the fixed-point tensors, 
and in an attempt to develop a theory of topological or- 
ders based on T- and S'-matrices, in this paper, we will 
calculate the T- and iS-matrices from the fixed-point ten- 
sors. This allows us to identify various topological orders 
in the the many-body wave functions described by the 
fixed-point tensors. In particular, we can calculate the 
number of quasiparticle types, the quasiparticle statis- 
tics and quantum dimensions, etc . 



II. LOCAL UNITARY TRANSFORMATION 
APPROACH AND MODULAR 
TRANSFORMATIONS 

In this section, we will review the local-unitary trans- 
formation approach and associated fixed-point tensor de- 



scription of many-body wave functions of topological or- 
dered states. The key in the local unitary transformation 
approach is a new definition of equivalence class on many- 
body ground states. Two gapped many-body ground 
states were defined as belonging to the same equivalence 
class {ie the same "phase") if and only if they are con- 
nected by a local unitary (LU) transformation defined as 
follows: 

|$(1)) ^ 1^(0)) iff 1^(1)) = r[e-'^o ^^^(^)]|$(0)) (1) 

where T is the path-ordering operator and H{g) = 
^iOi{g) is a sum of local Hermitian operators. It can 
be proved that the above definition is equivalent to the 
standard definition through "adiabatic evolution" {ie two 
states belong to the same phase if and only if they can be 
connected through a gapped adiabatic evolution).^ The 
new definition of equivalence class is advantageous be- 
cause it provides a very operational way to determine 
if two states belong to the same phase, thus provides a 
very natural way for renormalization. Under this defi- 
nition, "long-range entangled" states are those that are 
not in the same phase as direct-product states. Different 
"long-range entanglement" are then called "topological 
orders" . 

Following the new definition of equivalence class, 
a wave function renormalization scheme can be 
introduced. ^'^^'^^ It contains two basic moves called "F- 
move" and "P-move", both of which are the general- 
ized local unitary (gLU) transformations. Under this 
renormalization scheme, gapped ground states will only 
flow within the same equivalence class, ie renormalization 
does not change its topological order. Thus a fixed point 
in this renormalization scheme can be used to represent 
a whole equivalence class. Following the convention of 
Ref. 5, we represent the two renormalization moves at 
fixed point by tensors: 




also, the simplest fixed-point wave function was repre- 
sented by: 

*fix(^(7)) =^' = ^'*- (3) 

The above-mentioned fixed point tensors can be ob- 
tained through solving a set of self-consistent equations. 
Since the fixed-point state will not change under renor- 
malization, we can apply renormalization moves in ar- 
bitrary orders. This arbitrariness gives us a set of 
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self-consistent conditions that fixed point tensors must 
satisfy. The self-consistent conditions are summarized 
below: ^ 

Nikj = Nkji = Ni*j*k* > 0, ^^^Nii*kNjk*j* > 1, 
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(4) 



Conditions (4) form a set of non-linear equa- 
tions with variables A^ij/c, F^j^!^^^ ^ pkj,a(3 ^ 
{0f^0pi^0p2^0ai^0a2)^ which are tensor labels for fixed 
points. Following the previous discussion, the solutions 



(N^Jk. Fi^innx. AO win then give us a labeling of 

different topological orders. 

After obtaining such a tensor labeling for the ground 
state wave function, the next natural questions to ask 
are: Are they in one-to-one correspondence with different 
topological orders? How can we physically understand 
these tensors? Are there any physical quantities that we 
can at least numerically measure? 

The non-Abelian Berry's phases T and S are intro- 
duced precisely to answer these questions. Not only can 
T and 5* give a one-to-one label for different topolog- 
ical orders, but they can also provide a link between 
the fixed-point ground states and their corresponding 
quasi-particle excitations. We believe that the descrip- 
tion given by the T- and ^'-matrices is complete, meaning 
that they can completely characterize different non-chiral 
topological orders. Further numerical calculations of var- 
ious topological properties thus can all be obtained from 
T- and ^-matrices. 

Now we want to briefly introduce the concept of non- 
Abelian-Berry's-phase description of topological order 
which was first introduced in Ref. 8. The more detailed 
description will be given in the next section of the paper. 
The non-Abelian Berry's phases are obtained from two 
types of transformations named T- and S-transformation, 
both defined on a torus. The "T-transformation" can 
be defined by twisting a fixed-point graph along an axis 
(let us say the x-direction) by 360°; whereas the "S- 
transformation" can be defined by rotating a fixed-point 
graph by 90°. Essentially they are all discrete deforma- 
tions of a fixed-point graph on a torus. As hinted from 
the one-to-one correspondence between the ground state 
degeneracy on a torus and the number of different types 
of quasi-particles, the information of the quasi-particles 
is always encoded in the geometry of the torus. Here 
again, through the non-Abelian Berry's phases of the 
ground states, we can "decode" the information about 
their quasi-particles. So T- and ^'-matrices give a better 
description of topological orders. In this paper, we will 
label different topological orders by their resulting T and 
S'-matrices. 

Also mathematically, the relationship between T- and 
S'-matrix labels and fixed-point-tensor labels can be un- 
derstood within tensor category theory. On one hand, 
the T- and ^'-matrices correspond to quasiparticle exci- 
tations, which may have fractional and/or non-Abelian 
statistics and are described by a modular tensor category 
T. On the other hand, the fixed-point tensors correspond 
to the many-body wave functions of non-chiral topolog- 
ical orders, which are described by spherical fusion cat- 
egory J^. Thus to calculate the T- and iS-matrices from 
the fixed-point tensors is to calculate the quasiparticle 
excitations from the many-body wave functions, which 
corresponds to calculating the modular tensor category 
T from a spherical fusion category J^. Mathematically, 
calculating T from F corresponds to taking the "Drinfeld 
center" of the spherical fusion category Our calcula- 
tion of T- and ^-matrices from the fixed-point tensors 
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corresponds exactly to taking the "Drinfeld center" of 
the spherical fusion category T . 

The above mathematical picture agrees perfectly with 
our calculations: we perform modular transformations on 
several fixed-point states (III A, IIIB, IIIC, HIE, IIIF), 
all having the structure of spherical tensor category; our 
resulting T- and ^'-matrices for the quasi-particle exci- 
tations are all of the structure of modular tensor cat- 
egory. This result further strengthened our belief that 
tensor category theory is the mathematical structure be- 
hind topological order. 

In this paper, we present a lot more fixed-point solu- 
tions to eqn. (4), representing many more topological or- 
ders. We then apply "modular transformations" to many 
of the obtained results (III A, IIIB, IIIC, HIE, IIIF) to 
get their T and ^'-matrices, describing their correspond- 
ing quasi-particles. 

Note that there's a special solution, ie the "chiral" case 
III H. In that case, the labeling fixed-point tensors form 
the so-called "multi-fusion category" and its fusion rule 
breaks the "chiral" symmetry. Although the resulting 
fixed-point wave function has a trivial topological order 
without any symmetry, we anticipate that the "chiral" 
case can be highly non-trivial after introducing certain 
symmetries. 



III. LIST OF DIFFERENT TOPOLOGICAL 
ORDERS 

In this section, we present a list of all fixed-point 
solutions obtained so far from solving eqn. (4). The 
solutions consist of fixed-point gLU transformations 
(A/^j/e, P/^-^'"^) and fixed-point wave functions 

A\ For most of the cases (III A, IIIB, IIIC, HIE, IIIF), 
we also present their corresponding T and S'-matrices 
from modular transformations. 



A. iV = 1 loop states — the Z2 states 

Let us first consider a system where there are only 
two states |0) and |1) on each link of the graph. We 
choose = i and the simplest fusion rule that satisfies 
eqn. (4) (we call the N tensors satisfying eqn. (4) as "fu- 
sion rule"): 



A^ooo = A^iio = A^ioi = A^oii = 1, 
other Nijf. = 0. 



(5) 



Since Nijk < 1, there are no states on the vertices. So 
the indices a, /3, ... labeling the states on a vertex can be 
suppressed. 

The above fusion rule corresponds to the fusion rule 
of the N = 1 loop states discussed in Ref. 12, thus the 
name = 1 loop states. Note further that the three edge 
labels of N also form a Z2 group: for example, Nhq = 1 
represents the group action 1 1 = (Here our group 



action is addition). Thus we call the states obtained in 
this section the Z2 states. 

Due to relation eqn. (4), different components of ten- 
sor F^jl^ are not independent. There are only four in- 
dependent potentially non-zero components, denoted as 



r^OOO 
^000, 



^111, 



^011, 



77.110 
^110 



= /o 



/ 7^011 

i^ioo 



^001 



r = (j^oi°oM X)* 
= /i 

= /2 



(6) 



We note that F^jl^ in eqn. (2) relates wave functions 
on two graphs. In the above we have drawn the two 
related graphs right after the F tensor. The first graph 
following F corresponds to the graph on the left-hand 
side of eqn. (2) and the second one corresponds to the 
graph on the right-hand side of eqn. (2). The doted line 
corresponds to the |0) -state on the link and the solid 
line corresponds to the |l)-state on the link. There are 
four potentially non-zero components in P^^-^ , which are 
denoted by pov^Ps: 

Po^«=Po, Po''=Pu P?'=P2. P?'=Ps^ (7) 

We can adjust the total phases of Pi and to make 
Po > and > 0. We can also use a gauge transforma- 
tion to make /i > 0. 

The fixed-point conditions (4) form a set of non-linear 
equations on the ten variables /i, Pi, and A^. Many of the 
non-linear equations are dependent or even equivalent. 
Using a computer algebraic system, we simplify the set of 
non-linear equations. The simplified equations are (after 
making the phase choice described above) 

/o = /i = /2 = 1, fs r], 



Pi=P3 = Wo, P2 =Po, 

P0 + |Pl|' = l, |P2|' + |P3|' = 1, 

PrA'=P2A\ rjpsA'=prA', \A 



(8) 



0|2 



il|2 



1 



where rj = ±1. The above simplified equations can be 
solved exactly. We find two isolated solutions parame- 
terized by 77 = ±1: 

/o = /i = /2 = 1, fs Vi 



^0- ^ A' - 

V2' V2- 



(9) 



We also find 



e^^^2 = 1. (10) 
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Q- 



000 






FIG. 1: (Color online) The 4 characteristic non-contractable 
graphs on a torus. (The inner square and the outer square 
are the physical boundary of a torus.) Any other fixed-point 
graphs can always be reduced to the above four by F-moves 
and P-moves, but the four types cannot be changed into each 
other. The dotted line represents the |0)-state on the edges, 
and the solid line represents the |l)-state on the edges. The 
three integers ijk at the lower right corner of each graph are 
the edge labels following the convention of eqn. (11). 



The above solutions (F^]^;^^, i^^^'"^, A') with rj = ±1 
both correspond to fixed-point states. The 77 = 1 state 
is an equal weight superposition of all the graphic states 
that satisfy the fusion rule, whereas the r] = —1 state is 
also a superposition of all fusion-rule-satisfying graphic 
states, but some with coefficient 1 and some others with 
coefficient —1. 

Having obtained the abstract fixed-point solutions, we 
now want to physically understand the results by intro- 
ducing the concept of "modular transformations" . Mod- 
ular transformations are defined on a torus {ie a pla- 
nar graph with periodic boundary conditions in both x- 
and ^/-directions). Notice that after putting the above 
states onto a torus, there will be four different types 
of non-contractable graphs (see Fig. 1). They corre- 
spond to four types of fixed-point states that are lin- 
early independent (meaning they cannot be transformed 
into each other through F or P-moves) and form the 
degenerate ground-state subspace of a local Hamilto- 
nian (See Appendix A for the form of the Hamiltonian) . 
The modular transformations can be defined within this 
ground-state subspace. We will first introduce the "T- 
transformation" , also known as the method of Dehn 
twist. As in Ref. 14, we define the Dehn twist formally by 
requiring it to map a fixed-point state j, a, f3) = 



a/3 
fix 



/r 



on a torus to another fixed-point 



state defined on a different graph 4>^^(i, j, /c, a, / 
/I \\ 



a/3 
fix 



Then we can use an F-move to de- 



form the graph 




which leads to 



a unitary transformation T between the four fixed-point 
states on the torus, called the "T-transformation" . The 
deformation process could be seen more clearly from the 
following: 



a/3 
fix 



Dehn twist -f-a3 



1^5 




k J 



{F — move) 



{Deformation) 



(11) 



In the last step above, we have suppressed the vertex 
indices a, /3, (5, x because a maximum of one state is 
allowed on each vertex. 

Writing "T-transformation" in the basis of the four 
non-contractable fixed point states on the torus, we get 
a 4 by 4 unitary matrix with elements F^^j^^ . This matrix 
is called as "T- matrix". The basis of the 4 by 4 matrix 
is defined as follows: 



fix 




(1,0,0,0), 



fix 



(0,1,0,0), 



fix 



(0,0,1,0), 



*^>fix 




(0,0,0,1). 
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We now apply Dehn twist to all the four non- 
contractable states (recall Fig. 1) to get the T-matrix: 



fix 



000 



^000 ^fix 



fix 



^101 '^'fix 



Q- 



^fix 



I 



fix 



We thus obtain the 4 by 4 T-matrix for the Z2 states as 
follows: 











000 

Fio^i 

F^IS 

V F^ll 



(12) 



To complete the story, we will also introduce the "5- 
transformation" . Unlike T-transformation which is gen- 
erated by Dehn twist on the torus, iS-transformation 
is generated by 90° rotations. We define the 90° 
rotation formally by mapping a fixed-point state 



a/3 
fix 



on a torus to 



its clockwise 90° rotation state on a different graph 



a(3 
fix 



Then we can 



again use an F-move to deform the graph 



I a" 



back 



to 




, which leads to a unitary transformation S 



between the four fixed-point states on torus. This de- 
formation process can be seen more clearly through the 



following graph: 



ai3 
fix 



90 rotation 
>*fix 





{F — move) 



(13) 



Again, we have suppressed the vertex indices a, /3, ^, x 
in the last step above. 

Writing "^-transformation" in the same basis as "T- 
transformation" , we get a 4 by 4 unitary matrix with ele- 
ments F^i^j^ . This matrix is called as the "^-matrix" . We 

I J L 

now apply 90° rotations to all the four non-contractable 
states (recall Fig. 1) to get the iS-matrix: 



fix 



fix 



"101 



fix 



^fix 



77^000 /Fj 

^000 '^'fix 



i^/o°/*fix 



On 



I 



We can thus obtain the 4 by 4 ^-matrix as follows: 



F^ll 
Fl°l 
V OF/ 



(14) 



110 
10 



The T and ^-matrix obtained above together give us the 
modular transformations. 
When 7^ = 1, we have: 



a \ 

10 

1 

vo 1 y 



(15) 
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/I 

10 

10 

\0 1 



(16) 



We can choose a new basis to make the T and /S-matrices 
be in a more standard form. In this new basis, we require 
that the T-matrix to be diagonahzed and at the same 
time, ^-matrix to satisfy the fohowing requirements: 



1. M .-/c 



2' — ^Jiii 

3. Sii> 0. 



is a non-negative integer. 



(17) 



It is conjectured that the T and S'-matrices satisfying 
the above requirements are unique (up to permutations 
of the basis) , and thus the basis of T and iS-matrices are 
completely fixed through the above requirements (up to 
permutations.) Numerical simulations have been done 
in several cases in the later sections of the paper, all 
agreeing with the conjecture. 

In the above, Si^j represents the (i, j)-th matrix ele- 
ment of the ^-matrix, and S**^- is its complex conjugate. 
In Modular Tensor Category theory, these are precisely 
the conditions satisfied by the "modular matrix" S and 
the "diagonal twist" matrix T; thus the name for our T 
and S'-matrix. The first condition above is the so-called 
Verlinde Formula; it gives the fusion rule for the quasi- 
particle excitations. (Note that this is different from 
the fusion rule 5, which is the fusion rule for fixed-point 
ground states.) 

Now, diagonalize the T-matrix and make the iS-matrix 
satisfy 17, we have: 



/I 

10 

10 

VO -1 



(18) 



Note that the K = {^^ U{1) x U{1) Chern-Simons 

(CS) theory also has 4 types of quasiparticles with sta- 
tistical angles 

(6'"-) = (1,1,1,-1)- 

Therefore, the r] = 1 T2 state can be described by the 
above CS theory. 

When 7^ = — 1, we have: 



S: 



/I 








1 














1 


\0 


-1 





/I 











1 





1 








\0 


- 


-1 



(20) 



(21) 



Again, we can diagonalize the T-matrix and make the 
/S-matrix satisfy 17 under a proper basis; numerical sim- 
ulations show that such a basis is unique (up to permu- 
tations). In such a basis, we have: 



/I 

-i 

i 

Vo 1 



(22) 



S = 




(23) 



In this particular case, the above T and S matrices can 
be further reduced to the following forms: 



1 

i 



1 

-i 



(24) 




(19) 



Numerical simulations show that the above T and S- 
matrices are unique (up to permutation of basis.) The 
above T and ^-matrices match exactly with the corre- 
sponding matrices in the "Toric Code modular tensor 
category" in Ref. 15. From the dimension of the T and 
S*- matrices, we can see that there are four different types 
of quasi-particles. From calculating eqn. (17), we can 
get the fusion rules between the quasi-particle excita- 
tions. The quasi-particle statistical angles are given by 
the eigenvalues of T-matrix.^' The quantum dimensions 
of these four quasi-particles are given by the first row 
(column) of the ^-matrix. 



S 



^fl 1 
V2V -1 



^fl 1 
V2V -1 



(25) 



which shows the "doubled" structure of the quasi- 
particles. The above T and ^-matrices match exactly 
with the "doubled Semion" modular tensor category in 
Ref. 15. Again, we can see that there are four different 
types of quasi-particles from the dimension of the T and 
S'-matrices. From calculating eqn. (17), we can get the 
fusion rules between the quasi-particle excitations. The 
quasi-particle statistics are given by the eigenvalues of 
T-matrix, whereas the quantum dimensions of different 
quasi-particles are given by the elements on the first row 
of S'-matrix. 

Notice that the eigenvalues of T-matrix again match 
^0 2\ 

22] theory (which is equivalent 



with the K 
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to the K = theory), which has 4 types of 

quasiparticles with statistical angles 

Also, the "doubled" structure hints that the effective the- 
ory can be "decoupled", which is exactly what we get 

here: K = ("^ ^^j can indeed be decoupled. Therefore, 



the 77 = — 1 Z2 state can be described by the K 
CS theory. 



2 
2 2 



B. = 1 string-net state — the "Fibonacci" state 

To obtain another class of simple solutions, we modify 
the fusion rule to 

Nooo = Niio = Nioi = Noil = Nm = 1, 

other Nijk = 0. (26) 

while keeping everything the same. The above Nijj^ also 
satisfies eqn. (4). 

The new fusion rule corresponds to the fusion rule for 
the = 1 string- net state discussed in Ref. 12, thus 
the name = 1 string-net state. Notice the additional 
fusion rule can be written as 1^1 = 001 which looks 
like Fibonacci's golden rule, we also refer to the state as 
the "Fibonacci State". 

Again, due to the relation eqn. (4), the different com- 
ponents of the tensor Fj^j!^ are not independent. Now 
there are seven independent potentially non-zero compo- 
nents which are denoted as /o,...,/6: 



77^000 

^000 , 



77^000 

^111/ 

7^011^ 

^011/ 

77^011 

^111 



77^110 
^110 

77^110 
^111 

77^111 
^111 



= /o 






_ / 77^011 W 

-i^ioo,- 


Xr = 


/plOl 








= {fZH. 


Xr = 




= {FmH 


Xr = 





= {Fl^lH, X)* = h 



X=^Fl 



111 
10 



)* = /5 



(27) 



Note that F's described by fi and are the only F's 
that change the number of |l)-links and the number of 
|l)|l)|l)-vertices. So we can use the local unitary trans- 
formation Q^i^^ 1-^(1^ ^111) to make fi and /s to be positive 



real numbers. (Here Mi is the total number of |l)-links 
and Mill is the total number of |l)|l)|l)-vertices.) We 
also use the freedom of adjusting the total sign of F^jl^ 
to make Re(/o) > 0. 

There are five potentially non-zero components in P^^ 
which are denoted by pov7P4' 

p°^=po, p^'=pi, pr=P2, 

P?'=P3, Pl'=P4. (28) 

We use the freedom of adjusting the total phase of P^^ 
to make po to be a positive number. We can also use the 
freedom of adjusting the total phase of A* to make to 
be a positive number. 

The fixed-point conditions (eqn. (4)) form a set of non- 
linear equations on the variables /i, Pi, and A*, which can 
be simplified. The simplified equations have the following 
form 

^2 



/o — /l — /2 — /s — 1, /4 — /s — 



p\fl 



1, Po = hPi, P2 



/e >0, 

Po, P3=Pl, Pi 







iA'f = l, /| + /4-l = 0. 



(29) 

Let 7 be the positive solution of 7^ + 7 = 1: 7 = . 
We see that /s = y^. The above can be written as 



/i 

P2 



fo-- 
Po -- 

A° ^ 
We also find 



/2 — /s 

7 



7 



1, /4 = - 
Pl =P3 
1 



-/6 



7, 
1 



V7^ 



h = V7, 
= , P4 = 0, 

(30) 



The fixed-point state above corresponds to the = 1 
string-net condensed state^^ . 

Following the logic from last section, we will now ap- 
ply the modular transformations. Let us first consider 
Dehn twist on a torus. Note that we now formally have 
5 possible non-contractable graphs on a torus rather than 
4, due to the additional fusion rule. (Recall Fig. 1; the 
additional graph would be a graph with all solid lines.) 
Doing the Dehn twist for all 5 graphic states using (11), 
we have: 



fix 



^fix 



77^000 

^000 ^^'fix 



9 



fix 



fix 





101 / 












^1 



J^iio^fix 



O- 



We thus obtain the 5 by 5 T-matrix for the Fibonacci 
state as follows: 



000 



\ 

F^ll 

F^l^l 

F/io F.Yoi 



(32) 



V Fill / 

S 

/ 



Similarly, we can also get the S'-matrix by applying 90° 
rotations to all five non-contract able states on a torus: 




fix 




^fix 



^fix 



- ^000 '^'fix 



= ^^lVo°*fix 



o- 



cr. 



<0*fix 



From which we can get the 5 by 5 S'-matrix for the Fi- 
bonacci state: 



S ■ 



/ 7^000 





V 





7^110 

^110 




7^110 

^111 







z^lOl 
^101 








7^011 

^011 






\ 

77111 * 

^110 




FIR J 



(33) 



Now we can substitute in the values of the F tensors, and 
thus obtain: 



S = 



/I 























1 











1 











7 








V7 


Vo 


V7 








-7 


/I 

















7 








V7 











1 











1 








Vo 


V7 








-7 



\ 



(34) 



(35) 



where 7 



75-1 



There is, however, a small complication here compared 
to the previous Z2 case. Although there are now five 
possible non-contractable graphs on a torus, they are no 
longer linearly independent (meaning they can be trans- 
formed into each other through F or P-moves) and are 
not all fixed-point graphs. Thus they don't all corre- 
spond to ground states of a local Hamiltonian. As can 
be checked using the Hamiltonian construction in Ap- 
pendix A, the ground state subspace in this case is only 
4-fold degenerate. We thus need to restrict our modular 
transformations to be within this ground state subspace. 
In order to get the T and S'-matrix within the ground 
states, we need to project the above obtained matrices 
onto its 4 X 4 ground-state subspace. Such a process 
can be easily done by diagonalizing the Hamiltonian and 
projecting the T and ^'-matrices only to the ground-state 
subspace; see Appendix A for details. 

After the projection, we can again diagonalize the T- 
matrix and make the ^-matrix satisfy 17 under a proper 
basis; numerical simulations again show that such a basis 
is unique (up to permutations). We thus have the final 
form of T and 5*: 



/I \ 

e'^ 

e-'^ 

VO 1/ 



(36) 



V5 



10 



1 

2 

2 

3+75 
2 



1 + 75 
2 

-1 

3+75 
2 

-1-V^ 



l + ^/5 
2 

3+75 
2 

-1 

-1-75 



3+75 \ 
2 ^ 

-1-75 
2 

-1-75 
2 

1 



(37) 
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As can be easily seen, the above /S-matrix is real and 
symmetric. Also in this particular case, the above T and 
S matrices can be further reduced to the following forms: 

which shows the "doubled" structure of the quasi- 
particles. The above T and ^'-matrices match exactly 
with the doubled "Fibonacci" modular tensor category 
in Ref. 15. As before, we can tell there are four dif- 
ferent types of quasi-particles from the dimension of T 
and S] we can also get the fusion rules between quasi- 
particle excitations from calculating eqn. (17). The 
quasi-particle statistical angles are given by the eigenval- 
ues of T, whereas the quantum dimensions of different 
quasi-particles are given by the first row (column) of the 
5- matrix. 

Following the last section, here we want to also com- 
ment on the eigenvalues of T-matrix. These eigenvalues 
correspond exactly to the quasiparticle statistical angles 
of the doubled 5*0 (3) Chern-Simons gauge theory: 

(ei^') = (l,l,e^^,e-^^). 

Thus the "Fibonacci" state can be described by the dou- 
bled 5'0(3) CS theory. 



denoted as /o,...,/i3: 

77^000 . . . .. r 

^000. Jo 

7^000 ■•..../ '(^ / 77O22 W ( -^^202■.: \* 

=F^ < y:=h 

F^iiH X = iFzH Xr = h 
f'aX< X = {F',^X Xr = F^ X 
= iF???H Xr = u 

F^iiH X = (^^m H Xr = (F^o^X. Xr 

= Fi^^H X = h 
F^iiH X = (FiSiX. Xr = k 
F^ifX X = «2 H Xr = f^^Xk X 

= iF?^iX. Xr = h 

FllX i X = h 

FiiX i X=(^^noH Xr = h 



C. An iV = 2 string-net state - the "Pfaffian" state 

Here we will give a more complicated example of a non- 
orientable string-net state. We choose = 2, 0* = 0, 
1* = 1, 2* = 2, and 



A'ooo — A'oii — A/'iio 

= N112 = A^121 



A^ioi 

A^211 



^"022 
1. 



202 



^"220 
(40) 



-.112 
112 



nllO^ 
221 / 



[^121 
121 



n220 
220 



^-/lO 



>■■< X = (^2i^oH X)* = (^mH Xr 



X 



/ r.121 
1^210 

^111, 

/ 7^211 

(^211 



13 



X=/ii 

X)* = /: 



(41) 



The above Nijj. satisfies eqn. (4). 

Similar to the Z2 case, the three edge states of N form 
a structure called a "tensor category" with the follow- 
ing fusion rules between its three elements: I^cf = <j, 
l(g)^ = ^, ^(g)^ = l, cr(g)cr = 10?/^. In the above 
fusion rules, "1" represents state on the edge, "a" rep- 
resents state 1 and "7/^" represents state 2. (For example, 
N112 corresponds to a^a = ip.) Since the state also 
corresponds to the "Pfaffian" state in the Ising model, 
we have the name - "Pfaffian state" . 

Due to relation eqn. (4), different components of the 
tensor F^j^^ are not independent. There are fourteen 
independent potentially non-zero components which are 



There are ten potentially non-zero components in P^^ , 
which are denoted by pov^Po* 



pOO _ pui _ pU2 _ puu _ pui _ 



:>01 



j02 



.00 



:>01 



02 



:>00 



P7, 



:)01 



:)02 



= P9. 

(42) 



Using the "gauge fixing" discussed in section Ref. 5, we 
can consistently fix the phases of /i, /2, /4, /?, /g, /ii, Po 
and to make them positive. (Note that this 'gauge fix- 
ing' is non-trivial, since these phases are not completely 
independent; for e.g. the phases of /4, fj are complex 
conjugates of each other.) 



11 



The fixed-point conditions (eqn. (4)) form a set of non- 
linear equations on the variables /i, p^, and A\ which 
can be solved exactly. After applying the "gauge fixing" 
discussed above, we find only one solution 



/o — /i = 
/s = /g = 

Pl=P4= P8 



fr — fii — fi3 — 1, 
1 



no 



V2 



12 



-1, 



fix 



V2 



, P5 = 0, 



PO = P2 = P3 = P6 = P7 = P9 



2' 



We also find 



2' 72 



fix 



1. 



(43) 



(44) 



i^lVl^fix 



Now we are ready to apply the Modular transforma- 
tions. First we try to get the T-matrix by applying the 
Dehn twists. In the "Pfaffian" case we have altogether 
10 non-contract able graphs on a torus as can be seen in 
below. Applying Dehn twists to all of them using (11), 
we have: 




*^>fix 



^fix 



*^>fix 



*^>fix 



*^>fix 



^fix 



I 




^000 '^'fi^ 



Q- 




We thus obtain the 10 by 10 T-matrix for the Pfaffian 
state: 



/I 00000 00 0\ 

000010 00 

000000 00 1 

000100 00 

0^0000 00 ^0 

0^0000 00 

000000 -1 

000000 01 

00000 -1 00 

\0 01000 00 0/ 



(45) 



o- 



I 



Similarly, we can apply the 90° rotations and obtain the 
S'-matrix. Without going into to much details, we give 
here the final result (under the same basis): 



/I 00000 00 0\ 

0^0000 00 ^0 

001000 00 

000010 00 

000100 00 

000000 -1 

00000 -1 00 

000000 00 1 

0^0000 00 

\o 00000 01 oy 



(46) 




As in the "Fibonacci" case, the above 10 non- 
contractable graphs are not all fixed-point graphs and 
do not all correspond to ground states of a local Hamil- 
tonian. Using the Hamiltonian construction in Appendix 
A, it's easy to check that the true ground state subspace 



12 



is actually 9-fold; we can thus project the above matrices 
onto its ground state subspace. By carefully choosing the 
basis, we can again make the resulting 9x9 T- matrix di- 
agonalized and the 9x9 S'-matrix satisfy 17 at the same 
time, as shown below. Numerical simulations again show 
that such a basis is unique (up to permutations.) 



0* = 0, 1* = 1, 2* = 2, and 
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n 
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n n 




n 
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— if 
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n n 
u u 
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n n n 
u u u 
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— i 


n n 
u u 
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n n n 
u u u 
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U 


U ( 


3 8 U 
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U U U 
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(4 






































-10 





















e^'-f 




















lyi 






V2 


1 


V2 


2 


V2 


1 V2 


1 \ 


V2 





-V2 


2 





-2 


V2 




1 


-V2 


1 


V2 


-2 


V2 


1 


1 


V2 


2 


V2 











-V2 -2 


-V2 


2 





-2 











-2 


2 


V2 


-2 


V2 











2 


-V2 


1 


V2 


1 


-V2 


-2 


-72 


1 ^ 


1 


V2 





-V2 


-2 





2 


V2 


-V2 


Vi 


-V2 


1 


-V2 


2 




1 -^/2 


1 / 



(48) 

The above ^-matrix is again real and symmetric. Also 
in this particular case, the above T and 5* matrices can 
be further reduced to the following forms: 



(49) 















T= 




• 


^ 


e-*t 


VO 







VO 


-1 



1 



5=^1^/2 

2 M -V2 1 



y/2 1 
V2 



1 \/2 
•s/2 



000 



iVoii 
iViii 



A^iio 



A^ioi 



N21I 



N202 
1. 



N220 

(51) 



The above Nijk satisfies eqn. (4). 

These fusion rules correspond to the N = 2 string-net 
model with non-orientable strings in Ref. 12. Following 
that paper, we will also call this state the Ss state. As 
can be seen in the following, after gauge fixing, the fixed- 
point solutions completely agree with the "local rules" 
obtained in Ref. 12, thus we re-obtain another string-net 
state from local unitary transformation point of view. 
(Note that edge labels 1 and 2 are reversed in Ref. 12 
compared to this paper.) 

Due to relation eqn. (4), different components of the 
tensor F^j!^ are not independent. There are fourteen 
independent potentially non-zero components which are 
denoted as /o,...,/i8: 



77^000 

^000 , 



ZT^OOO 

^111/ 



ZT^OOO- 
^222 ;• 



^011. 



;^011 



^122 



i.=fo 



■k - 1^200 .■ 



\* / z?202 
.) -1^020, 



7^220 
^002 



(50) 



{F'2^H. Xr 

{F^iri Xr 



77^111 
^011 



77^121 

^012 



J* 



which shows the "doubled" structure of the quasi- 
particles. The above T and 5*- matrices match exactly 
with the doubled "Ising" modular tensor category in 
Ref. 15. As before, we can tell that there are 9 types 
of different quasi-particles from the dimension of T and 
5* and through calculating the Verlinde Formula 17, we 
can get the fusion rule between them. The quasi-particle 
statistics can be obtained from the eigenvalues of T- 
matrix and the quantum dimensions of these quasiparti- 
cles are given by the elements on the fist row of /S-matrix. 



^211 



7^022 \ 
^022 , 



^111 



- ^011 
X=(^2\^ 



202 
202 / 



/ 7^211 

1^102 



- i^lOl 
= /6 
*=/7 

'Xr = F^ii 

*=/8 



X)* 



D. Another N — 2 string-net state — the state 

By adding an additional fusion rule on top of the pre- 
vious result, we can get yet another example of a non- 
orientable string- net state. Here we still choose N = 2^ 



^110 



77^110 

^111 

7^110 

^112 



111 

10 



112 

10, 



J* = /lO 

J* = hi 



13 



p220\. 
-^^220/ 



/l2 
/l4 

(^^2^1^: 

{Fill 

fl8 



Xr = /i3 

Xr = Fi?iHX 
Xr = /i5 



J* = {F?20 



X = he 

Xr = hr 



(52) 



There are ten potentially non-zero components in P^^-^ , 
which are denoted by pov^PQ* 



dOO 


Dll 

:)02 



P5, 



Pi 



■Pi, 
-P6, 



:>02 



p02 



P2, 
P7, 



00 



pOO 



P3, 



01 



pOl 
^2 



P2'^=P10. 



= P4, 

(53) 



Using the "gauge fixing" discussed in section Ref. 5, we 
can fix the phases of /i, /2, /s, /s, /lo, /ii, /i6, Po and 
to make them positive. Again, these phases are not 
completely independent, so the "gauge fixing" has to be 
done self-consistently. 

The fixed-point conditions (eqn. (4)) form a set of non- 
linear equations on the variables /i, p^, and yl*, which 
can be solved exactly. After applying the "gauge fixing" 
discussed above, we find only one solution 



/o — /l — ••• — /s — /l6 — /l7 — /l8 — I5 
1 

2' 

1 

2 

1 

P0=P2=P3=P7 =P8= PlO 



/g — fii — /i4 — 

flO = -fl3 = 
Pl=P4= P9 



fl2 = 0, 
/l5 = -1, 
P5=P6= 0, 



We also find 



1 

7H' 



2 

7i* 



(54) 



(55) 



A careful comparison with Ref. 12 shows that the F- 
tensors obtained above perfectly match the F local rules 
for the N = 2 non-orientable string-net state in Ref. 12 
(after switching the edge labels 1 and 2). Thus the above 
state corresponds to the standard lattice gauge theory 
with gauge group Ss - the permutation group of 3 ob- 
jects. 



E. An N — 2 string-net state — the Z3 state 

The above four examples all correspond to non- 
orientable string- net states. Here we will give an exam- 
ple of an orientable string- net state. We choose N = 2, 
0* = 0, 1* = 2, 2* = 1, and 

A^OOO = A^012 = A^120 = A^201 = A^021 = A^102 = N210 

= A^iii = A^222 = 1. (56) 

The above Nijk satisfies eqn. (4). 

This state also appeared in Ref. 12 as the other ex- 
ample of N = 2 string- net state. As in the Z2 case, the 
three edge labels of N form a Z3 group after we switch 
the definition of positive direction in the third label of N. 
For example, N012 becomes A^oii after switching, which 
will then correspond to group action 0(^ 1 = 1. we will 
call such a state the Z3 state. 

Due to relation eqn. (4), different components of the 
tensor F^j!^ are not independent. There are eight inde- 
pendent potentially non-zero components which are de- 
noted as /ov J7: 



7^000 .. . . 

^000/ 




= /o 






F??rr{ 


X 




Xr 


= F^S^ < X. 






= 


X)* 


= fi 


E^ooo'-..../ 

^222 / \ 


C 


= (F^SS"^., 


T)* 


= iF^?^r\ X^ 








X= 


/2 




X 


= F^ilH 


X= 


(^^2^o°2^H Xr 






= (fXK 


Xr 


= h 




X 


= (Fl^lH 


Xr 


= F^ X 






= (Fl^iH 


Xr 


= h 




X 


= (Fi^iH 


Xr 


= ^01^2)^ X 






= (FiSlH 


Xr 


= h 


FSsH 




= {FM?)i 


Xr 


= (Fii^H X) 



^120 rx ^ 



Fli^}{ 



iF 



210\. 



220/ 



= /6 



(57) 



There are nine potentially non-zero components in P^^-^ , 
which are denoted by pov^Ps* 

/'0°°=P0, Po°'=Pl, Po°'=P2, Pr=P3, P{'^=P4, 
Pr=P5, P2°°=P6, P'2'=V7, P^'=PS. (58) 

Using the "gauge fixing" discussed in section Ref. 5, we 
can fix the phases of /i, /2, /e, Po and to make them 
positive. 
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The fixed-point conditions (eqn. (4)) form a set of non- 
linear equations on the variables /i, p^, and A\ which 
can be solved exactly. After applying the "gauge fixing" 
discussed above, we find only one solution 

/, = 1, i = 0,l,...,7, 
1 



S'-matrix under a change of basis: 



Pi i 0,1,..., 8, 

73 



(59) 



We also find 



(60) 



The fixed-point state corresponds to the Z3 string-net 
condensed state^^. 

The process of applying the Modular transformations 
are exactly the same as in III A. By applying Dehn twists 
to all the 9 non-contract able fixed-point states on a torus, 
we can get the T-matrix. Note that as in III A, here all 
the 9 non-contractable graphs correspond to fixed-point 
states and are ground states of a certain local Hamilto- 
nian. Without going into too much details, we present 
here the resulting 9 by 9 T-matrix: 

/1 00000000\ 
010000000 
001000000 
000001000 
T = 000100000 (61) 
000010000 
000000001 
000000100 

Vo 00000010/ 

Similarly, S'-matrix can be obtained by applying the 90° 
rotations to all the 9 fixed-point states. The resulting 9 
by 9 ^-matrix is then: 
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0000 














0001 


[0 








1 


0000/ 



As before, we can choose a particular basis in which 
T-matrix is diagonalized and S'-matrix satisfies eqn. 17 
in modular tensor category theory. Without going into 
too much details, we present here the resulting T and 
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(63) 



(64) 



1 \ 

r 
1 

(0' / 



where ^ = e*^. As before, all information of the quasi- 
particles can be obtained from the above T and S matri- 
ces. 

Note that the corresponding eigenvalues of T are: 
(l,l,l,l,l,e'^,e'^,e-^*,e-'^). 

which exactly correspond to the statistical angles of the 
U{1) X U{1) Chern-Simons theory^^.i^ 



(65) 



with K = ^3 • This is also equivalent to Z3 gauge 
theory. 

F. N = 3 string-net states — the Z4 states 

Now we will allow four states on each edge of the graph, 
namely |0), |1), |2) and |3). We will first give an example 
of orientable string-net state. We choose TV = 3, 0* = 0, 



N220 
N211 
(66) 

The above Nijk satisfies eqn. (4). 

We will see that one of the solutions to be obtained 
corresponds to the Z4 gauge theory. Again, as in the 
Z3 case, the three edge labels of N form a Z4 group 
after switching the positive direction of the third label of 
N: for example, N121 becomes N123 which corresponds 





1* = 3, 2* = 2, 


3* = 1, 


and 






(62) 


A^ooo = A^oi3 


= ^'l30 


= A^301 


= Nq22 


= N202 = 




= A^osi 


= A^103 


= A^310 


= Nii2 


= N12I = 




= N233 


= A^323 


= Nss2 


= 1. 
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to group action 102 = 3. Thus, we will call states 
obtained from the above fusion rules the Z4 states. 



Again, due to relation eqn. (4), different components of 
tensor F^jl^ are not independent. There are now twenty 
independent potentially non-zero components which are 
denoted as /o,...,/i9: 



nOOO 

000 , 



000 
11 . 



222 , 



r.000 
333 . 



r^Oll 
Oil 



fo 

/pOllV^ 

l-f'300, , 



-^003/ 



(^o3or\ = /i 



1-^200. 
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77^220 
^002 
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E^033 
^322 



X 



7^112 
^132 



7^112 
^310 



7^123 
^210 



E.130 
^110 



n220 
220 



X 



fiflM X : 



7T.231 

^013 



K X=iFi^iH Xr 



(F, 



i321 
301 



Xr 
(Fi^iH Xr 
(Fio'^H Xr 
/12 

(^^3^l2°)-( Xr 



fio 

Fi^^H X 
hi 

(fh^X Xr 



E^332 
^130 



H 



/: 



13 



(FiiXi Xr 



7T1231 

^120 



i310 
223 



Xr = /i4 



'^'H Xr 



1^213 

FfuHX 



(F 



231 , 



715 



(Fi^SH 

(Fi^iH 

{Fii^H 

fl8 
fl9 



Tie 



ri220 

113 



717 



(X 



(67) 



Oil 



F 



7^011 

^233 



E^022^ 
^022. 



E.022 
^133 



r^022 
^311 



^ - (^2^03^; 

= iF?^i 

= ^^0^3\^K 

; = (FisiX. 

X = (-^232 : 

X = (Fio'i 
= FiiH 



Xr = F, 
Xr = h 
Xr = (F, 
X=/< 



123 
032 



123 
103 



HXr 



* = /7 

Xr = F^iiy, 
Xr = fs 
Xr = (F 
X=/! 



^.303 
212 



H Xr 



There are now sixteen potentially non-zero components 
in P-^-^ , which are denoted by pov^Pis- 

pOO _ ^ pOl _ ^ p02 _ ^ p03 _ ^ pOO _ ^ 
— PO, ^0 — PI5 ^0 — P2, ^0 — P3, ^1 — P4, 

pOl p02 p03 _ pOO pOl _ 



:)02 



p02 
^3 



PlO, 
^Pl4, 



:)03 



p03 
^3 



^Pll, 
^Pl5- 



:)00 



^Pl2, 



Pl3, 



(68) 

Using the "gauge fixing" discussed in section Ref. 5, we 
can fix the phases of /i, /2, fs, /s, /s, /g, /ii, /i3, /i4, 
/i7, and to make them positive. Again, these 
phases are not completely independent, so the "gauge 
fixing" has to be done self-consistently. 

The fixed-point conditions (eqn. (4)) form a set of non- 
linear equations on the variables /i, p^, and A\ which 
can be solved exactly. After applying the "gauge fixing" 
discussed above, we find two solutions 

/12 = /l5 = /l9 = ^, 
/i = 1, i = otherwise, 
1 



Pi = 2' ^ = 0, 1,...,15, 



(69) 
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where rj = ±1. We also find 



we have the resulting T and S- matrices: 



(70) 



/ 1 



Similar to the Z2 and Z3 case, we can obtain the T 
and S'-matrices of the two states above using modular 
transformations. By relating the eigenvalues of the T- 
matrices to quasi-particle statistics, we can also find the 
corresponding effective theory for these two states. First, 
by applying Dehn twists to all the 16 fixed-point states 
on a torus, we obtain the following T-matrix: 



/I 



(71) 



V I 



Similarly, we can apply the 90° rotations and obtain the 
16 X 16 S'-matrix on the same basis: 



T = 



(73) 



-1/ 



/I 1 
1 1 

1111 
1111 
1 1 
1 1 
1 1 
1 1 
1 -1 
1 -1 
1 -1 

% 1 

\-\ i -i \ 
\-\ i -i -1 
\-\ i -i -1 
Vl -\ i -i \ 



llllllllllllllX 
111111 -1-1-1-1-1-1-1-1 
1 —1 —1 —1 —i—i—i—ii i i i 
1 —1 —1 —1 i i i i —i —i —i —i 
1-11-1-11 1-1-11 1-1-11 
1 —1 —1 1 1 —1 —ii i —i i —i —i i 
1 —1 —1 1 1 —1 i —i —i i —i i i —i 
1-11-1-11-11 1-1-11 1-1 
i i 1 —i i —1 1 —i i —1 1 i —i —1 



1 i —i 1 —i —1 1 i i —1 1 —i 
1 i —i 1 i 1 —1 —i —i 1 —1 i 
i i —1 —1 i —i 1 —1 —ill 
—i —1 1 i —i —1 1 —i i —1 
i i 1 i —1 1 —i —i —1 1 i 
i i 1 —i 1 —1 ill —1 —i 
—i —1 —1 —ill —1 i —i \ I 

(74) 



/I 



s = 



(72) 



All information of the quasi-particles can be obtained 
from the above T and S matrices. 

Note that the eigenvalues of T-matrix are: 



(1, 1, 1, 1, 1, 1, 1, 1, -1, -1, -1, -l,i,z, -i, -i) 



which correspond exactly to the statistical angles of the 
Z4 gauge theory, or equivalent ly the /7(1) x /7(1) Chern- 
Simons gauge theory: 



(75) 



As before, we can choose a basis so that T-matrix is 
diagonalized and S'-matrix satisfies eqn. 17. For 77 = 1, 



with K 



4 

4 oy 



On the other hand, when 77 = —1, by properly choosing 



17 



the basis, we have the following T and S- matrices: 

/ 1 \ 



T = 



e 4 



(76) 



1 / 



y^llllllllllllllllX 

11-1-11 1-1-11 1-1-11 1-1-1 

1 —1 i —i 1 —1 i —i 1 —1 i —i 1 —1 i —i 

I —1 —ill —1 —ill —1 —ill —1 —i i 

II 1 1 1 1 1 1-1-1-1-1-1-1-1-1 
11-1-11 1-1-1-1-11 1-1-11 1 
1 —1 i —i 1 —1 i —i —1 1 —i i —1 1 —i i 

1 1 —1 —ill —1 —i i —1 1 i —i —1 1 i —i 
4 1 1 1 1 —1 —1 —1 —1 —i—i—i—ii i i i 
1 1 —1 —1 —1 —1 1 1 —i —i i i i i —i —i 
1 —1 i —i —1 1 —i i —i i 1 —1 i —i —1 1 
1 —1 —i i —1 1 i —i —i i —1 1 i —i 1 —1 
1111 —1 —1 —1 —1 i i i i —i —i —i —i 
1 1 —1 —1 —1 —1 1 1 i i —i —i —i —i i i 
1 —1 i —i —1 1 —i i i —i —1 1 —i i 1 —1 
\1 —1 —i i —1 1 i —i i —i 1 —1 —i i —1 1 / 

(77) 

In this particular case, again, the above T and S matrices 
can be further reduced to the fohowing forms: 



/I \ 

0-10 

e^^ 

Vo 



4 



/I 

0-1 
e-'- 
\0 



(78) 



or equivalently the U{1) x U (1) Chern-Simons gauge the- 
ory: 



£ = 1-K"ai^d,ajxe''''\ 



(81) 



with K = 



4^ 
4 4/- 



G. More N — 3 string-net states — the Z2 x Z2 states 

Now we will give an example of a non-orient able string- 
net state with four edge states. That is to set = 3, 
0* = 0, 1* = 1, 2* = 2, 3* = 3. We then have the 
following solution for N: 



000 



= A^oii 


= A^ioi 


= A^iio 


— -^022 


— N202 — 


= -^033 


= -^303 


= -^330 




= N231 = 


= iVi32 


= A^321 


= ^"213 


= 1. 





N220 

N312 

(82) 



The above Niji^ satisfies eqn. (4). 

As before, we can recognize from the above that the 
three edge labels of N form a Z2 x Z2 group if we flip 
the positive direction of the third label. Since the state is 
non-orient able, switching the direction has no real effects. 
Note that any element square will give unity, which is 
precisely the property for group Z2 x Z2. 

Similarly, due to relation eqn. (4), different compo- 
nents of the tensor F^jl^ are not independent. There are 
now twenty-two independent potentially non-zero com- 
ponents which are denoted as /ov5/2i- 




(79) 



which shows the "doubled" structure of the quasi- 
particles. As before, all information of the quasi-particles 
can be obtained from the above T and S matrices. 
Note that the eigenvalues of T now become: 



(1, 1,1,1,1,1,-1,-1, e±'4,e±^T) 

which correspond exactly to the statistical angles of the 
U{1) X U{1) Chern-Simons gauge theory 



4iv\ 



(80) 



r.OOO 
^000 . 



^111.. 



000 



r^OOO 
^333 , 



7^011 ^ 

^011, 



/o 








/ 77^011 W 

i^ioo, 


Xr = 


i^oioXX ^ 


^)* 


^001/ \ 


y=h 






/ r.022 W 
1^200, 


Xr = 


/ 7;;.202 

(^020 XA ^ 




7^220\.../ 
^002 / ■.. 


y=h 






/ E^033 W 
1^300, 


Xr = 


/ r.303 

1^030 ^ / 




7^330 \ . . . 
^003/ •■. 










Xr = 


U 
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in P-^-^ , which are denoted by Pov?Pi5* 



77^011^ 

^233 . 



77^011^ 

^322, 



E^022^ 
^022, 



F 



O22W 

133 . \ 



r^022\ 
^311, 



77^033 
^033 , 



^122, 



^211/ 



7^110^ 

^110/ 



77^110^ 

^223/ 



(F^^f < Xr = h 
{FisiX. Xr = iF?^i 



77^321 

^012 



SX=/. 



/ 77202 

1^202 / 



*=/7 



- ^02^3^ ^ - (^302^) 



( z?303 
1^212 



101 

232 



77312\ 



( Z?303 
1^303 



77^123 

^032 



/ 77.312 
1^203 



101 
323 



103 



13 



*=/8 

* = {F'.ll'r 
= h 



Jw 



{FmH Xr 



y 

Xr = hi 



iF^?!H 



77^213 

^031 



..)* = /: 



12 



3* = (F?^^) 



J* 







■PO, P^'=PU P^'=P2, P^'=PS, Pi 
Pr=P5, P?^=P6, Pr=P7, P^°=P8, 
P^'=P10, P^'=PU, P^'=P12, Pt=PlS 



>00 



pOl 
^2 



:)02 



Pl4, 



^03 



Pa, 



(84) 



Using the "gauge fixing" discussed in section Ref. 5, we 
can fix the phases of /i, /2, /s, /e, /s, /g, /ii, /12, 
/i47 /167 Po and to make them positive. As before, 
these phases are not independent to each other, so the 
gauge fixing has to be done in a self-consistent way. 

The fixed-point conditions (eqn. (4)) form a set of non- 
hnear equation on the variables /i, p^, and A% which 
can be solved exactly. After applying the "gauge fixing" 
discussed above, we find four sets of different solutions, 
listing below: 

/. /, = !, i = 0,l,...,21, 
1 



Pi 



^0 = ^1 



0,1,..., 15, 



IF fl3 — /l5 — fl7 — /18 — /l9 — /20 

/i = 1, i = otherwise^ 



A° = -A' 



-A' = A-" 



-1, 



PO = P3 = P4 = P7 = P8 = Pll = P12 = Pl5 



Pl=P2=P5=P6=P9= PW = Pl3 = Pl4 



2' 



2' 



332/ 



FlMr 

^220/ 



7^220 
^113 



132 
310 



(Fi 

77^330 
^112 



F^iiH X = (F^li 



118 



X = (FiliX 



^r = (F, 

C = /l5 

X)* = /16 

= /l7 



312 
130 



in. /i3 = /i5 = /i7 = /21 = -1, 

= 1, i = otherwise^ 

1 

2' 



Po=P2= Pa 

Pl=P3= P5 



^ PlA 



^ Pl5 



1 

"2' 



4O 



-A' 



IV. fis = /l9 = /20 = /2I = —1, 

/i = 1, i = otherwise^ 



^331/ 



(F 



231 
320 



77^330 
^221 



FiilH X = (Fiil 



(Fii^H 



119 



Xr = /20 



J* 



P0=Pl=Pi=P5=PS=P9= P12 = Pl3 = 
P2=P3=P6=P7= PlO = Pll = Pl4 = Pl5 

AO = = -A^ = -A^ = -. 



2' 



1 

' 2' 
(85) 



77330\. 
^330/ 



/21 



(83) 



We also find 



There are now sixteen potentially non-zero components 



(86) 
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We can further check if the results obtained above are 
consistent with the ones in the Z2 case (c.f. Ill A). A 
quick check can be done as fohows: In III A, we have 

/o = 1 /3 = T]. 

(We use f to differentiate it from the / obtained in this 
section.) Since Z2 x Z2 can be seen as two copies of Z2, 
we expect 

/o = fS' X fo = 1, /i3 = /a' X f^' = v\ 

/18 = fS' X /I' = r?^ /21 = fr X /I' = r?" X r?^ 

Thus /i3, /18, /21 would either all be one, or two out 
of them would be —1. This precisely corresponds to the 
four solutions obtained above. A more detailed check 
also shows consistency between the two results. 



H. An N — 3 string-net state — the "Chiral" state 

Finally, we present here a new state which was once 
thought to be a non-string-net state (but eventually 
turned out not to be one, fortunately or unfortunately). 
In this state, interestingly enough, the chiral symmetry 
is spontaneously broken. As before, we choose = 3, 
0* = 0, 1* = 3, 2* = 2, 3* = 1, and 

A^ooo = A^013 = A^130 = A^SOl = ^"123 = A^231 = A^312 

= A^222 = 1. (87) 

The above Nijk satisfies eqn. (4). 

This state has several interesting features. Besides its 
chiral symmetry being spontaneously broken, the three 
edge labels of N form a structure called a "multi-fusion 
category" . This mathematical structure has the following 
fusion rule between its four elements: 

Mij X Mki = SjkMii, where ij,k,l = 1 or 2. 

After making the following mapping, Mn = state 0, 
M12 = state 1, M22 = state 2, M21 = state 3, we can 
see that the fusion rule above is exactly equivalent to 
(87) (again, we need to flip the positive direction of the 
third label of N). Although the fusion rule is fairly non- 
trivial, we will see that at least in the no-symmetry case 
which we are studying in this paper, the corresponding 
fixed-point state is trivial. 

Due to relation eqn. (4), different components of the 
tensor F^j!^ are not independent. There are now six inde- 
pendent potentially non-zero components which are de- 



noted as /ov J5: 



z^OOO ■ , . . ■ _ 


fn 
J u 






^ 333,- \ ..-V 


V-^ 100 . 




(Fn?n j^ : V 
V-' 010 / A j,-^/ 




7^310\../ 
-''^001/ 


y=/i 




^ 233,- \ ..-V 


323 / \ 




-'011/ \ 




/ 7^323\^ 


J.) - 


J2 


^no^ ^ = 


(^3^32°) ( 


Xr- 


h 








F???HX 




{FiiiH 


Xr = 





FiiiH X = /5 (88) 

There are eight potentially non-zero components in P^^-^ , 
which are denoted by pov^Pr* 

pOO _ ^ p03 _ ^ pOO _ ^ p03 _ ^ p21 _ ^ 

Fq —Po, Fq —Pi, —P2, ^1 —Ps, ^2 —Pa, 

P2^=P5, Pt=P6, PF=P7- (89) 

After a gauge fixing process,^ we can fix the phases of /i, 
/3, /4, Po and to make them positive. 

The fixed-point conditions (eqn. (4)) form a set of non- 
linear equation on the variables /i, Pi, and A*, which can 
be solved exactly. After applying the "gauge fixing" dis- 
cussed above, we find the solution with an undetermined 
parameter rj: 

/, = !, z = 0,l,...,5, 

P0=P2=PA=P6= 

pi=P3=P5=P7 = \/r^-^, 

A^ =A^ = r]^l-r]^. (90) 

We also find 

It can be shown that the above state corresponds to 
a trivial loop state with adjustable loop weights. By in- 
troducing the "double-line rule" and associate different 
weights to the "dashed-line" loops & "solid-line" loops, 
the state reduced to a state with independent loops. This 
can be seen very clearly in Fig. 2. 

IV. SUMMARY 

In this paper, we computed the T and ^'-matrices that 
were introduced to define topological orders in 2D.^'^ We 
argued that these matrices could be used to fully describe 



20 



State 
state 1 
state 2 
state 3 



Fusion rule: 

N222 




□ 




FIG. 3: (Color online) Both graphs represent non- 
contract able graphs on a torus. The graph on the right is 
a zoom-in version of the graph on the left. When construct- 
ing the Hamiltonian, we use the graph on the right to show 
the details of F and P- moves. 



FIG. 2: The "double-line rule" . Replace the single-line triva- 
lence graph by double lines. If different link states are defined 
as shown above, then the fusion rules are automatically satis- 
fied. (The fusion rule N222 = 1 was shown in the graph as an 
example.) Suppose we associate a weight 77 to a dashed-line 
loop and a weight ^1 — rf to a solid-line loop respectively, 
then the weight of any graph can be determined. (Note that 
all the loops are independent, making it a trivial loop state.) 



topological orders. We also showed explicitly how these 
matrices could be obtained from their fixed-point states 
through "modular transformations". The final matrices 
were uniquely determined after a proper choice of basis. 

More specifically, The "modular transformations" were 
defined through two transformations, ie the "Dehn twist" 
and "90° rotation" . After applying both transformations 
in the ground- state subspace, we got two matrices, T 
and describing "modular transformations". We then 
changed the basis to get the proper form of these two 
matrices. 

To apply the two transformations in the ground-state 
subspace, we used projections. We first defined T and S 
for all non-contractable graphs on a torus, then applied 
unitary transformation to diagonalize the Hamiltonian. 
We then projected T and S only onto the subspace with 
the smallest eigenvalue of the Hamiltonian, thus obtained 
T and S for the ground states. 

To get the "proper form" of the two labeling matri- 
ces, we first diagonalized T- matrix, and then searched 
for unitary transformation that didn't change the diag- 
onal form of T, and at the same time, made S satisfy 
a few requirement including the "Verlinde" formula (see 
eqn. (17)). After this unitary transformation, the form 
of T and S were uniquely determined only up to permu- 
tations of basis. 

The resulting T and ^'-matrices we obtained in this pa- 
per are unique labels for different topological orders. We 
believe they contain full information of the phase and 
can fully characterize different topological orders. Fu- 
ture work can be done to further generalize this method 
to obtain T and S matrices for general tensor-product 
states or other many-body states, so as to determine the 



topological order of an arbitrary many-body state. Even 
further, using this T and ^'-matrices description, we may 
eventually be able to describe phase transitions between 
different topological orders. 

We would like to thank Z.-X. Liu, L. Kong for some 
very helpful discussions. This research is supported by 
NSF Grant No. DMR-1005541, NSFC 11074140, and 
NSFC 11274192. Research at Perimeter Institute is sup- 
ported by the Government of Canada through Industry 
Canada and by the Province of Ontario through the Min- 
istry of Research. 



Appendix A: Ideal Hamiltonian for fixed-point 
states on a torus 

The Hamiltonian construction is similar to the one in 
appendix of Ref. 14. The basic idea is to construct a 
Hamiltonian which is a sum of commuting projectors and 
is thus exactly solvable. On specific lattices for example 
the honeycomb lattice, it can be shown that all fixed- 
point wave functions {Nijk, F'^^i;"^:;^/ , P^^'''^ , A') that we 
obtained from solving eqn. (4) are exact gapped ground 
states of such a local Hamiltonian.^^ 

The Hamiltonian is of the following form: 

F = ^(l-Q„) + ^(l-Bp) (Al) 

V p 

where sums over all vertices and sums over all 

plaquettes. The Hamiltonian H acts on the Hilbert space 
Vq formed by all the graph states. Operator in H acts 
on the states of the 3 links that connect to the vertex v: 




Clearly, Q^j is a projector = Qv Operator Bp in H 
acts on the states of all the links and vertices belonging 
to the same plaquette p. 
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For our purpose in this paper, we only need to con- 
sider the Hamiltonian acting on the subspace of ah non- 
contractable graphs on a torus (See Fig. 3). The Bp 



operator on such a torus can be constructed from a com- 
bination of the F- moves and P- moves as fohows: 
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t,s i' k' 
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t,s i' k' j'* 
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* it* rpi' s* i rpk'* s ' 
ssi' ^ kj*k' ij' 



k* T^j s*j 
^i'k'*i' 




( TDts\* T^i*it* T^i' s*i T^k'* s*k* T^j' s* j j^isi' ^ 
2^ 2-^ 2-^ 2-^ 2-^\^i ) ^ssi' ^kj*k'-^j*ij'* ^ i'k'*i"-^i"*si"'*^ 
t,s i' k' j'* i" i'" 
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k\ 













\^\^\^\^\^\^\^fTDts\^T^i*it*T^i's*iT^k'*s*k*-pj's*j j^isi' Z?*'"*s*i*^ 
2^ 2-^ 2-^ 2-^ 2-^ 2-^ 2-^y^i ) ^ssi' ^kj*k'-^j*ij'* ^ek'*e'-^i"*si'"*-^k'*j'k"*^ 
k' j'* i" i'" k" 



t,s 
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t,s i' k' j'* i" i'" k" j" 



J Si 



TTil S 



s z pk'* s*k* rpj • 



J 771252 7712 

i"-^i"*si"'*-^k' 



''*s*i* rpk" s* k' ^ 





\s \ 1 
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t,s i' k' j'* i" i'" k" j" 




_ \ ^ \ ^ \ ^ \ ^ \ ^ \ ^ \ ^ \ ^( TDts\^ T^i*it* T^i's*i T^k'*s*k* i^j' s*j j^isi' T^i"*s*i* j^k" s* k' ^ 
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/ 



_ \ ^ \ ^ \ ^ \ ^ \ ^ \ ^ \ ^ \ ^ \ ^( -nts\* T^i*it* T^i's*i T^k'*s*k* T^j's*j j^isi' T^i'"*s*i* j^k" s*k' ^ 

- 2^ 2^ 2^ 2^ 2^ 2^ 2^ 2^ 2^y^i ) ^sse ^kj-k'-^j-ij'- ^i'k'H"-^i"''si"'--^k'-j'k"--^j'i"*j"^ 

t,s i' k' j'* i" %"' k" j" t' 




thus we have the matrix B given by 

^ - /^z^z^z^z^z^z^z^z^y^i ) ^ssi' ^kj^k'^jHj'^ ^i'k'-i"^i"^si"'^^k'^j'k"^^j'i"^j"^ 

t,s i' k' j'* i" %'" k" j" t' 

Ff:g;:,Fi;,f;/" p^y (A4) 



which is the matrix form of operator Bp. 

Recah that ah fixed-point states are exact ground 
states of the Hamiltonian H defined in eqn. (A2). Thus 
by finding the exact ground states of H among ah the 
non-contractable graph states, we can find the subspace 
of ah the fixed-point states. 

In sections IIIB and IIIC, we are required to reduce 
the modular transformations to only within the ground- 



state subspace of the above Hamiltonian. This can now 
be easily achieved. We can first define T and 5* on all 
the non-contractable graphs on a torus, then find the 
ground-state subspace by diagonalizing the Hamiltonian 
on these non-contractable graphic states. We can then 
project T and S only onto the ground-state subspace of 
the Hamiltonian, thus obtain the required T and S. 
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